Influence of solvent quality on polymer solutions: a Monte Carlo study of bulk and 

interfacial properties 
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The effect of solvent quality on dilute and semi-dilute regimes of polymers in solution is studied by 
means of Monte Carlo simulations. The equation of state, adsorptions near a hard wall, wall-polymer 
surface tension and effective depletion potentials are all calculated as a function of concentration 
and solvent quality. We find important differences between polymers in good and theta solvents. In 
the dilute regime, the physical properties for polymers in a theta solvent closely resemble those of 
ideal polymers. In the semi-dilute regime, however, significant differences are found. 
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I. INTRODUCTION 

The properties of polymers in solution are determined 
by the balance between effective monomer-monomer at- 
tractions and excluded volume repulsions. Upon increas- 
ing the relative strength of the attractions, for example 
by cooling, several regimes are encountered. As long as 
the thermal energy fc^T far exceeds the absolute value e 
of the effective monomer- monomer attraction, good sol- 
vent conditions prevail, and individual polymer coils are 
swollen due to the dominance of excluded volume effects. 
The polymer radius of gyration Rg scales like L'^ , where 
L is the number of monomers or (Kuhn) segments, and 
V K, 0.59 is the Flory exponenti*2,. When the tempera- 
ture is lowered, the coils shrink due to the action of effec- 
tive (solvent-induced) monomer-monomer attraction, un- 
til the theta regime is reached; at the theta-temperature 
Tg, monomer-monomer repulsion and attraction cancel, 
at least at the two-body level, and individual coils behave 
essentially like ideal (random walk) polymers, such that 
their radius of gyration scales like . Below Tg individ- 
ual coils collapse into dense globules with Rg ^ Ls . 

Moving away from the infinite dilution limit, inter- 
actions between different polymer coils come into play, 
and give rise to the lowest order correction to the van 
t'Hoff limit of the osmotic equation of state, valid for 
non-interacting polymers (see Eq. ^ below). The cor- 
rection, due to pair interactions, is proportional to the 
square of the overall monomer concentration c; the coef- 
ficient is the second virial coefficient B2{T; L). The Boyle 
temperature Tb is the temperature at which B2 vanishes, 
i.e. B2{Tb,L) = 0. Since this condition is met when the 
monomer-monomer repulsion and attraction cancel, one 
expects that Tb — Tg. 

Finally, upon lowering the temperature at a finite poly- 
mer concentration, the polymer solution is found to sep- 
arate into polymer-poor and polymer-rich phases below 
a critical temperature Tc. Within Flory-Huggins mean- 
field theory Tc is found to coincide with Tb within cor- 
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FIG. 1: Schematic phase-diagram of a polymer solution as 
a function of temperature T and monomer concentration c, 
for a finite length L. Rg decreases with decreasing T so that 
for good solvents, the monomer concentration at which the 
dilute regime crosses over to the semi-dilute regime increases. 
As T is lowered further, there is also a crossover to the "theta 
regime", denoted schematically by the shaded region. In the 
limit L ^ 00, Tb and Tc both approach Te, so that the size 
of the "theta regime" decreases. Furthermore, the monomer 
concentration at the crossover between the dilute and semi- 
dilute regimes, as well as that of the critical point, both tend 
to zero for infinite L. 



rections of order \/\/L. Thus, a polymer solution is 
characterized by three temperatures Tg, Tb, and Tc, of 
which the first is a single polymer property, defined in 
the L 00 limit. For any given L, it is believed that 
Tb(L) >Tg > Tc{L)^'^, but in the scaling fimit L 00, 
Tb = Tc = Tg. In fact, when simulating simple models of 
non-ideal polymers, probably the most accurate estimate 
for Tg in the scaling limit is obtained by extrapolating 
results for Tb{L) to 1/L — (fi>. In this paper the oper- 
ational definition of Tg will be Tg — Tb for sufficiently 
long chains {L > 500). 

For a finite L, there is a "theta solvent regime" of 
temperatures around Tg, where the polymer behavior 
most resembles that predicted for theta polymers. The 
longer the polymer, the sharper the transition between 
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"theta solvent" and "good solvent" regimesii^i^. A typ- 
ical phase-diagram for solutions of polymers of finite 
length is presented in Fig. ^ 

In the low concentration limit, the physical proper- 
ties of interacting polymers in solution closely resemble 
those of ideal polymers with the same Rg. For good sol- 
vents, however, important deviations from ideal behavior 
rapidly set in for increasing concentration, even in the di- 
lute regime, see e.g. refi^. On the other hand, the fact 
that B2 (T, L) ~ for T near Tg suggests that in the theta 
regime the osmotic pressure 11 will follow van t'Hoff's law 



(3U^^ (1) 

for a larger concentration range than is found for poly- 
mers in good solvent. Nevertheless, as concentration is 
increased further, the effect of higher order virial coeffi- 
cients will eventually kick in, leading to deviations from 
the simple van t'Hoff behavior. Similar differences with 
ideal-polymer behavior should also be observable for in- 
terfacial properties such as the adsorption or surface ten- 
sion near a hard non-adsorbing wall. 

While much theoretical and numerical work on simple 
models of interacting polymer solutions has been devoted 
to behavior for athermal conditions, 01 at T — Tg, less 
is known about how properties of polymer solutions vary 
with temperature and concentration in the intermediate 
regime (00 > T > Tg). 

In the present paper we present MC data for the bulk 
and interfacial properties of dilute and semi-dilute poly- 
mer solutions for the same model as used in refe, over a 
range of temperatures between the athermal (high tem- 
perature) limit and Tg. The main objective is a system- 
atic investigation of polymer solution properties over a 
wide range of concentrations, to determine how equilib- 
rium properties like the osmotic pressure, the adsorp- 
tion near a wall, the surface tension and the polymer 
induced depletion potential between two walls vary as 
thermal conditions change from the good solvent to the 
theta-solvent regime. In particular we are examining the 
important problem of how the behavior of polymers in 
theta-solvent differs from that of ideal (non-interacting) 
polymers as their concentration increases. 

The paper is organized as follows. The polymer lat- 
tice model and the computational methodology are in- 
troduced in section II. Results for the osmotic equation 
of state and semi-dilute regimes will be presented and 
discussed in section III. We shall next turn our attention 
to monomer density profiles near a non adsorbing wall as 
well as the polymer- wall surface surface tension as a func- 
tion of temperature and concentration (section IV). The 
polymer-induced depletion potential between two walls 
will be discussed in section VI, and concluding remarks 
will be made in section VII. 



II. MODEL AND SIMULATION 
METHODOLOGY 

A. Lattice model of polymer solutions 

A familiar coarse-grained representation of linear poly- 
mers in a solvent, which captures the essential physical 
features, is a self-avoiding walk (SAW) lattice polymer, 
where each lattice site can be occupied by at most one 
monomer (or polymer segment) to account for excluded 
volume, and where pairs of non-sequential monomers 
of nearest-neighbor (n.n.) sites experience an attrac- 
tive energy —f3e. The parameter e accounts for the dif- 
ference between solvent-solvent, solvent-monomer, and 
monomer-monomer interactions, and includes both ener- 
getic and entropic componentsi*^. As the dimensionless 
ratio pe = e/ksT increases from zero (corresponding to 
the athermal solvent limit) to larger values, the quality 
of the solvent decreases, i.e. effective monomer-monomer 
attractions become increasingly important, so that the 
polymer coils tend to shrink until the theta-regime is 
reached. Upon further cooling at finite concentrations, 
phase separation sets in at some T < Tg, see e.g. Fig.^ 

In this paper we consider polymer chains made up of L 
monomers in a simple cubic lattice of M sites (coordina- 
tion number z = 6), with periodic boundary conditions. 
If N polymers live on that lattice, the polymer concen- 
tration is p = N/M, while the monomer concentration is 
c = NL/M. Rg ^ L", (with v 0.59 in good solvent, 
and = i in theta solvent) is the radius of gyration of 
a polymer coil, and the overlap concentration is defined 
as 1/p* = ^nRg, where Rg, which depends on f3e, will 
conventionally be chosen to be the radius of gyration in 
the infinite dilution limit (p 0) limit, separates the 
dilute regime of the polymer solution {p/ p* < 1) from 
the semi-dilute regime {p/p* > 1). 

In the scaling limit (L 00), the properties of a poly- 
mer solution in the dilute and semi-dilute regimes de- 
pend only on p/p* and Rg, and are independent of the 
monomer concentration oi*S. In other words, simulations 
with different c but the same p/p* should give equiva- 
lent results when length is expressed in units of Rg. This 
requires that simulations be carried out with sufRciently 
long polymers so that the monomer concentration c <C 1. 
Under good solvent conditions, Rg « 0.4L°'^^ for poly- 
mers on a simple cubic lattice, so that 
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c* = Lp* « (2) 
whereas in theta-solvent, Rg ~ O.SSLs'^:, and so 

Most of the subsequent MC simulations have been car- 
ried out for L = 500 polymers, such that c* ~ 0.027 
under good solvent conditions, while c* ~ 0.06 in theta- 
solvent. In order to ensure that c < 0.25 in all simu- 
lations, the semi-dilute regime which could be explored 



was hence restricted to p/p* < 10 in good solvent, and 
p/ < A in theta-solvent. For higher concentrations, the 
c-dependence of the results would no longer be negligi- 
ble, and non-universal effects would become important, 
depending on the property under consideration. For this 
reason, much smaller polymers, say L = 100, could not 
be used to study the semi-dilute regime. 

The previous extensive MC study by Grassberger and 
Hegger— determined Tb as a function of L for a model 
identical to the one used here. They found that for L = 
500, the Boyle temperature Tg gave /3e — 0.265, the 
value which we shall use throughout the discussion in this 
paper as a reasonable estimate for the theta-temperature. 
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B. Monte Carlo simulation methodology 

We used 4 different types of moves to sample config- 
uration space: The pivot algorithm attempts to rotate 
part of the polymer around a randomly chosen segment. 
Translations move the whole polymer through space, and 
reptations remove a segment at one end of the polymer 
and attempt to regrow it at the other end. Since trans- 
lation and pivot rapidly become less efficient for increas- 
ing polymer concentration, we also used configurational 
bias Monte Carlo (CBMC) moves''. These are also very 
effective at lower temperatures, since the attraction be- 
tween polymers reduces the accepted number of pivot 
and translation movcs^. 

Simulations were carried out in a box of size 160 x 100 x 
100 for polymers of length L = 500. For calculations 
of radii of gyration, and other bulk properties, periodic 
boundary conditions were used in all directions, whereas 
for surface properties calculations, hard walls were placed 
at and 160 along the x-axis. The osmotic pressures were 
calculated using the method of Dickman^. A external 
repulsive potential that acts between on monomers 
next to the wall is introduced, and varied such that the 
parameter < A = e~^"/^'^'^ < 1. A=0 prevents any 
particles from being adjacent to the wall, and the varying 
wall hardness allows the the osmotic pressure /3n to be 
calculated from the integral 

/3n(p) = -pM), (4) 

where Pu,(A) is the site occupation fraction of monomers 
right next to the wall. We performed 5 simulations at dif- 
ferent A values, chosen from a standard Gauss-Legendre 
distribution, to carry out the integral above. By using 2 
walls of varying hardness, the statistics are enhanced over 
using a single wall. Adding repulsive walls increases the 
bulk density at the center of the simulation cell, an effect 
which becomes more pronounced for smaller boxes. Thus 
N/V should not be taken as the bulk density. Instead, 
we have determined the bulk density by averaging over 
the region, near the center of the box, where the influence 
of the two depletion layers have completely decayecilS^. 



FIG. 2: Effective exponents z/, extracted from simulations 
with different lengths L, as a function of fie. . Note that 
the value presented at /3e = 0.256, corresponding to Tb for 
L = 500, was calculated at slightly different temperatures, 
corresponding to Tb, for each polymer length. 



As a first application of the MC code, we have com- 
puted the radius of gyration of an isolated polymer for 
lengths L — 100 to L = 5000, and extracted an expo- 
nent Rg ^ L'^ , shown in Fig.|21 For higher temperatures, 
1/ = 0.588 within error bars, as expected for the good 
solvent regime, but for /3e = 0.2 the exponent v appears 
to deviate slightly, and we find v = 0.570 ± 0.003. In the 
limit i — > oo, this scaling should be independent of tem- 
perature, as long as T > Tg. However, here we observe 
some rounding, most likely due to finite length crossover 
effects. The precise temperature at which theta-type 
scaling behavior of Rg sets in is dependent on length. 
The longer the polymer, the smaller the range of tem- 
peratures over which transition from the good-solvent to 
the theta regime occurs^iii. Note also that for a given L, 
the Rg for e.g. /3e = 0.1, is smaller than Rg for f3e — 0. 
Nevertheless, the scaling with length, Rg ^ L'^-^^, is still 
the same for those two temperatures, i.e. the polymers 
are still in the good solvent regime. Thus the effect of 
lowering the temperature is simply to renormalize the 
effective step-length. 

For T = Tb, with Tb for each length taken from refj^, 
we find ly — 0.498 ± 0.008, as expected for the theta 
regime^^. As the temperature is lowered further, f con- 
tinues to decrease, and for low enough T the polymer 
should collapse to a compact globule state where ^ = \- 
The beginning of this trend is evident in the plot for 
fie — 0.285. The larger error bars reflect the more pro- 
nounced finite L effects expected at this temperature. 
For short polymers the transition from extended to col- 
lapsed states occurs over a broad range of temperatures, 
whereas for longer polymers the collapse transition is 
sharper^. 

Rg shrinks (for fixed L) as the temperature is low- 
ered, because the excluded volume interactions are par- 
tially compensated by attractive interactions between 
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with concentration for polymers in good solvent, is virtually 
independent of density for T = Tb, and increases with con- 
centration when T < Tb- For each temperature simulations 
were performed up to a maximum monomer concentration 
c ~ 0.25, so that the total p/ p* range accessible decreases 
with temperature. 



monomers. Similarly in the good solvent regime, Rg de- 
creases with density because the excluded volume inter- 
actions are screened by overlapping polymer coils. This 
decrease follows a scaling law Rg ^ psiS once the 
semi-dilute regime is well developed^^. At the theta- 
temperature, on the other hand, chain statistics are 
nearly ideal on all length scales and concentrations^,, so 
that Rg should be independent of concentration. For 
even lower temperatures, Tb > T > Tc, the screening 
of attractive interactions now implies that Rg should in- 
crease with concentration. These trends are indeed ob- 
served in our MC simulations, as illustrated in Fig.Ofor 
L — 500 polymers. 



III. EQUATION OF STATE IN THE DILUTE 
AND SEMI-DILUTE REGIMES 



FIG. 4: MC data for the reduced internal energy U(p)/e = 
—Nc for L = 500 chains, plotted as a function of polymer 
concentration p/ p* for 5 different temperatures. Note how 
the curvature changes with solvent quality. 
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FIG. 5: Linear graph of the EOS Z = l3Yl[p)/p as a function 
of concentration for different temperatures. (Inset) A blowup 
of the dilute regime highlights the differences between near 
theta solvents and polymers in good solvent. The former show 
behavior that mimics that of ideal polymers. 



The thermodynamic property which is most readily 
extracted from the simulations is the internal energy of 
polymer solutions, which is simply U = —Net, where 
is the average number of nearest-neighbor "contacts" be- 
tween non-connected monomers. The energy per polymer 
is U /N = —eLhc, where fic is the average number of non- 
connected nearest neighbors around a monomer. Within 
a simple mean-field picture, fic is expected to increase 
linearly with the polymer concentration. Such linear be- 
havior is indeed observed in Fig. 0] for theta conditions. 
At higher temperatures, there is a significant downward 
curvature on U{p), while below Tb, but above Tc, U{p) 
curves upward. 

The osmotic equation of state (EOS) Z = 
Ii{p) / {pksT) was calculated using Dickman's method^ 
described in the previous section. The MC results for Z 



as a function of p/ p*, for various temperatures, are plot- 
ted in Figures and El on linear and logarithmic scales. 
As expected, Z increases with p/ p*, and decreases with 
temperature^^. Since the theta-solvent regime is defined 
by T « Tb, the EOS is expected to be very flat at low 
polymer concentrations since Tb is defined as the tem- 
perature at which the second virial coefficient B2{T,L) 
vanishes in the virial expansioni^: 

Z=^ = l + B2C + B3c''+Oic^). (5) 
P 

This is indeed confirmed by the MC results in the di- 
lute regime, as illustrated in the inset of Fig. El Up to 
p/p* ~ 1, hardly rises with p, i.e. the system behaves 
like a solution of ideal polymers. At higher concentra- 
tions, however, interactions between monomer triplets 
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FIG. 6: The EOS as a function of concentration for differ- 
ent temperatures, plotted on a double logarithmic scale. In 
the semi-dilute regime, the EOS is consistent with scaling 
theory, which predicts /3n(p)/p ~ p^ "^ for good solvents and 
(3Yl{p)/ p ~ p^ for theta solvents. 

come into play, so that, according to Eq. jSJ, ^ is ex- 
pected to increase like p^, in agreement with a standard 
scaling argumentiiS. The same scaling argument pre- 
dicts that the EOS of polymers in good solvent scales 
as Z ~ pi/(3='-i) the exponent 1.3 indeed agrees 

with the slopes of the double logarithmic plots of Z{p) 
shown in Fig. H for p/p* > 1 and < /3e < 0.2. The 
slope of the theta-temperature EOS is clearly signifi- 
cantly larger, and compatible with the expected expo- 
nent of 2. However, for the L = 500 chains used in 
the present simulations, the accessible semi-dilute regime 
of concentrations is too small to allow a fully quantita- 
tive estimate (c.f. the discussion in section II). Since the 
EOS of theta-polymers scales with a larger exponent than 
that of polymers in good solvent, there will eventually 
be a concentration p/p* where the EOS of the theta- 
polymers becomes larger than that of polymers in good 
solvent. This may, at first sight, seem counter-intuitive, 
but it should be remembered that for a given size 
theta-polymers are more compact. By extrapolating our 
simulation data to larger p/p*, we expect the cross-over 
between theta-polymers and /?e = (SAW) polymers to 
occur at p/p* ^ 35 which is, in practice, very hard to 
achieve both in simulations and experiments^^. On the 
other hand, if one starts from a solution under good sol- 
vent conditions, i.e. "high" temperature, at a given p, 
and lowers the temperature, then the EOS will decrease 
monotonically, because Rg, and with it the ratio p/p*, 
decreases with temperature. 



IV. MONOMER DENSITY PROFILES NEAR A 
NON-ADSORBING WALL 

For entropic reasons, polymer solutions will be highly 
depleted in the vicinity of a non-adsorbing hard wall: 



since fewer conformations are available to a polymer coil 
the wall will effectively repel the polymers so that their 
concentration will drop from the bulk value for distances 
z > Rg, to a much lower value at contact-. The case of 
polymers in good solvent (in the SAW limit) has been 
investigated in much detail in ref.®, where important de- 
viations from ideal polymer behavior were found even in 
the dilute regime. Here we continue these investigations, 
but for polymers in solvents of varying quality. 



A. Monomer density profiles 

Fig [3 compares the present MC results for the re- 
duced monomer density profiles c{z)/c as functions of 
the distance z from the planar surface in the good sol- 
vent (/?£ — 0) and theta-solvent (/3e — 0.265) limits for 
L = 500 chains, and several bulk concentrations. The dif- 
ferences between the two regimes are very significant. In 
the dilute regime the density profiles in the theta solvent 
hardly change, and are remarkably close to that of ideal 
polymers, while the profiles in good solvent already differ 
significantly from the latter, and are "pushed" closer to- 
ward the hard wall. In the semi-dilute regime, on the 
other hand, the density profiles of polymers in theta- 
solvent differ considerably from that of ideal polymers 
(which is independent of concentration) and instead more 
closely resemble those of polymers in good solvent, al- 
though the depletion "hole" is systematically wider for 
the former compared to that of the latter. This is il- 
lustrated in Fig. ISl where the profiles at three different 
temperatures are compared at low concentration, and in 
the semi-dilute regime. 



B. Reduced adsorption 

The reduced adsorption of polymers near a hard wall 
is defined asi&: 

P Jg 

where 51^^ /A is the surface excess grand potential per 
unit area A, p is the chemical potential, and h{z) — 
c{z)/c— 1. Because of the depletion hole, r(p) is neg- 
ative, and has the dimension of length. One could re- 
place c{z) by the distribution of the center of mass (CM) 
of the polymers pcm{z)- Although the reduced profiles 
hcm{z) = Pcm{z)/p— 1 would look different from h[z) 
(see e.g. Fig.l of Refi^), the reduced adsorptions would 
be identical to those calculated from the monomer den- 
sity profiles, due to the conservation of the number of 
monomers. 

For ideal polymers, the reduced adsorption can be ex- 
actly calculated to bei2*ii: 

t"^ = -2Rg/V^ w -1.1284i?g, (7) 
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FIG. 7: Monomer density profiles c(z)/c for polymers near a 
hard non-adsorbing wall. The upper panel shows the density 
profiles for polymers in good solvent, while the lower panel is 
for polymers at T = Tb- At these near theta conditions, the 
density profiles in the dilute regime are nearly indistinguish- 
able from those of ideal polymers, but important deviations 
emerge in the semi-dilute regime. 



which is independent of concentration. For polymers in 
good solvent, a renormalization group (RG)^- calcula- 
tion predicts that r(0) ~ — 1.074i?g for the adsorption 
in the low-density limi^, which is close to, but slightly 
less negative than the ideal result As the polymer 
concentration increases, interactions between monomers 
of non-ideal polymers cause the depletion layer to shrink, 
as illustrated by the density profiles shown in Figs. [7| and 
IHl Thus, in contrast to the case of ideal polymers, where 
c(z)/c, and therefore F is independent of concentration, 
for interacting polymers the reduced adsorption F should 
become less negative for increasing polymer concentra- 
tion. This is indeed what is observed in Fig.|51 where the 
variation of the reduced adsorption with concentration is 
shown for several temperatures. 



FIG. 8: Monomer density profiles as a function of temper- 
ature for two densities. The upper panel is for p/p* « 0.14, 
and the lower panel is for p/p* « 4, well into the semi-dilute 
regime. The density profiles at lower temperatures are more 
extended when plotted in terms of zjRg. 
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FIG. 9: Reduced adsorption V jRg as a function of concen- 
tration for several temperatures. (Inset) Blowup of the dilute 
regime. Note the different curvature at low densities. The val- 
ues at p/p* = are take from theory to be F/iJg = 1.074— for 
the first three temperatures, and T / Rg = 2/-^?? for T = Tb- 



1. Adsorption in the dilute regime 

In the dilute regime {p/p* ^ 1), there is, once again, a 
qualitative difference between polymers in theta or good 
solvent, as highlighted in the inset of Fig. O While the 
curvature of the F(p) curves with concentration is nega- 
tive in good solvent conditions, it is clearly positive under 
theta-conditions, and changes to negative in the semi- 
dilute regime. Up to p/p* — 1, the results for F remain 
close to (but above) the ideal polymer result iQ. Note, 
however, that the relative deviation from ideal polymers 
is more pronounced than for the osmotic EOS (compare 
the insets to Figs. |Sland|5Jl, hinting at a stronger effect 
of three-particle interactions on the adsorption of theta- 
polymers. 



2. Adsorption in the semi-dilute regime 

In the semi-dilute regime we expect the adsorption 
to be proportional to the correlation length £,/Rg ^ 
(p/p*)"''/'-^'^"^-!^, since this is the only relevant length- 
scale. Note that the numerical prefactors for different 
semi-dilute correlation lengths vary with the property 
one is attempting to describe, see e.g42i for a very use- 
ful discussion of this matter. In our previous papes^ 
we showed that for polymers in good solvent, the re- 
duced adsorption T/Rg ~ {p/ p*)~°''^'^ is consistent with 
this scaling. Here we present more simulations in the 
semi-dilute regime to confirm this behavior, finding that 
T/Rg ~ — 0.4(p/p*)~° '''' provides a good fit. Scahng the- 
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ory considerations also predict that ^/i?g ~ {p/p*)~^ for 
theta solventsi*^. However, this behavior cannot be reH- 
ably extracted from the current simulations of F, which 
are not performed for a large enough range of p/ p*. Fur- 
thermore, we expect that for the higher values of p/ p* 
sampled here, finite c effects may already be coming into 
play. 

V. WALL-POLYMER SURFACE TENSION 

The wall-polymer surface tension, i.e. the free-energy 
cost of introducing a non-adsorbing hard wall and its 
associated depletion layer, can be calculated from the 
adsorption and EOS by use of the Gibbs adsorption 
equationSS: 

By performing one integration by parts w.r.t. density, 
Eq. ||HJ| can also be expressed as: 

7.(p) = -Jl{p)np) + (^^^ dp'. (9) 

The first term in this equation has an appealing physi- 
cal interpretation as the free energy cost per unit area of 
creating a slab cavity of width F(p). For ideal polymers, 
where F(p) is independent of concentration, this term is 
the only one that contributes, and so 7^ — —pT^'^ — 
— p2/y^. Since the EOS and T{p) for interacting poly- 
mers in various quality solvents were calculated in the 
previous sections, we can now, using Eq. determine 
the surface tension for different solvent qualities. Re- 
sults are shown in Fig.^] For good solvent conditions, 
i.e. fie — 0, they are, to within simulation errors, in 
near quantitative agreement with the RG calculations of 
Maasen, Eisenriegler, and Bringer^^. In fact, the larger 
number of simulation results in the semi-dilute regime, 
used in the present work, lead to even better agreement 
than that shown earlier in Ref.^. We are not aware of 
theories of comparable quality for the surface tension of 
theta polymers. 

As discussed earlier for the EOS and the adsorption, we 
again find qualitatively different behavior for the T = Tb 
solution in the dilute regime, as highlighted in the inset 
of Fig.^l The surface tension for such theta like solvents 
resembles that of ideal polymers much more closely than 
that of polymers in good solvent. As expected from the 
related behavior of the adsorption, the relative deviation 
from ideal polymer behavior in the dilute regime is some- 
what larger than what was found for the EOS. Note that 
the curves ^w/lid cross at very low p/ p*. This follows 
from the fact that the low concentration limit of the sur- 
face tension is given by limp^o 1^1 w = — pF, and — F is 
smaller for interacting polymers than for ideal or theta 
polymers. (See Fig. |5Jl. At slightly higher densities, the 
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p/p 

FIG. 10: Wall-polymer surface tension divided by the ideal 
surface tension for various solvent qualities, plotted as a func- 
tion of concentration p/p*. (Inset) A blowup of the dilute 
regime highlights the differences between ')whm for good and 
theta solvents. 

curves cross because 7^, increases more rapidly with con- 
centration for good solvents than for theta solvents. 

In the semi-dilute regime, where F ^ ^ ~ p-''/(3i'-i) 
and n(p) - p3../(3<.-i)^ p simplifies^ to 

7^^(p)«^nf^p2./(3.-i)_ (10) 

This implies 7^'* ~ pi-539 semi-dilute scaling for poly- 
mers in good solvent, and ^ f? behavior for theta 
solvents, which is consistent with the results in Fig. EH 
Furthermore, this scaling suggests that for large enough 
p/p*, the surface tension curve for theta polymers will 
eventually cross that of polymers in good solvent. In- 
deed, Fig. ^1 already appears to suggest this behavior, 
although some care must be taken in extracting the p/p* 
at which crossover occurs, due to the the possibility of 
finite c effects for theta polymerai^. With these caveats 
in mind, our results suggest that the surface tension for 
theta polymers would be greater than that of interacting 
polymers when p/p* > 10, which is considerably lower 
than the predicted crossing for the EOS, and within reach 
of experiment and simulations. On the other hand, just 
as in the case of the EOS, if one were to lower the tem- 
perature of an experimental polymer solution at a given 
p, this would lead to a lower p/p*, so that a monotonic 
behavior of 7u, with temperature is expected when fol- 
lowing this route. 

VI. DEPLETION POTENTIALS BETWEEN 
TWO WALLS 

The depletion potential between two walls or plates 
induced by a polymer solution is defined as the difference 
in free energy between the cases where the plates are at a 
distance z, or infinitely far apart. The polymer solution 
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between the plates is taken to be in equilibrium with 
a much larger reservoir of pure polymer solution at the 
same chemical potential. At infinite distance, the free 
energy cost of having two plates is simply twice the cost 
of making a depletion layer, while at contact, these two 
depletion layers are destroyed. Thus, when z — 0, the 
free energy per unit area of plate is given by 



W{0) = -27^(p). 



(11) 



For small z virtually no polymer is expected to pen- 
etrate between the plates, which are hence pressed to- 
gether by the osmotic pressure n(p), leading to a lin- 
ear increase of W{z) with z. As discussed in reS^, the 
simplest approximation for the depletion potential is to 
continue this linear form for larger z: 



W{z) = W{Q) + n{p)z- 
W{z) = ; 

where the range is given by 

W{0) 



z < D jp) 

z > D„{p) (12) 



n(p) n(p) ■ 



(13) 



We note that this approximation is similar to that 
adopted by Joanny, Leibler, and de Gennes who, in their 
pioneering paper^24, approximated the force between two 
plates as constant for x < 7r^(p) and zero for x > 7r^(p). 
This also results in a linear depletion potential. 

In ref^^, we showed that this simple theory is virtu- 
ally quantitative for polymers in good solvent^. Fur- 
thermore, it is well known to be quite accurate for ideal 
polymers as welU^. Although one could easily generalize 
the theory to include the small curvature seen for ideal 
polymers close to z = Dw^'^ , we ignore these small cor- 
rections in the interests of simplicity. Since this theory is 
accurate for polymers in good solvent as well as for ideal 
polymers, we postulate that the same simple assumptions 
are valid for other solvent qualities, and use Eq. H12|l to 
calculate the depletion potentials. 

The well-depth W{0) follows from Eq. (^TJ, and is 
shown in Fig. ^] The behavior mirrors that of Fig. ^| 
of course. From the inset it is clear that theta polymers 
most closely resemble ideal polymers in the dilute regime, 
while important deviations are found in the semi-dilute 
regime. 

In Fig. US we plot the range -Dm(p) for a number of 
different solvent qualities. In the dilute regime, high- 
lighted in the inset, the range for the near theta solvent 
is fairly close to that of ideal polymers, whereas for poly- 
mers in good solvent the range decreases markedly in the 
dilute regime. In the semi-dilute regime both good and 
theta solvent regimes show important deviations from 
ideal polymer behavior. From scaling theory we expect 



-0.77 



and Dii;(p) p for the polymers in 



good and theta solvents respectively. 

The results for W{Q) and Dm{p) can be combined with 
Eq. H12|) to determine the full depletion potential between 




FIG. 11: The depletion potential at contact is given by 
iy(0) = —2jw{p)- (Inset) A blowup of the dilute regime 
shows that theta polymers most closely resemble the ideal 
polymers. 
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FIG. 12: Range of the depletion potential Dw{p), as given 
by Eq. I|13^ . In the dilute regime, highlighted in the inset, 
theta polymers are most similar to ideal polymers, whereas 
polymers in good solvent already show fairly strong devia- 
tions: at p/p* = 1 the range has dropped by almost a factor 
2 compared to ideal polymers. 



two plates. An example, for p/p* w 1, is shown in Fig. 1131 
The depletion potential for theta polymers most closely 
resembles that of ideal polymers. For lower p/p* we ex- 
pect an even closer correspondence. 

In the dilute regime. Figs. I11H13I suggest that the de- 
pletion potential between two plates, induced by theta 
polymers, resembles that of ideal polymers. We have 
recently shown how to construct depletion potentials be- 
tween two spheres from that between two plates^^, and 
found quantitative agreement with direct simulations of 
L — 500 SAW polymers between spheres, and also good 
results for ideal polymers. This suggests that the same 
procedure should work well for the depletion potential 
between two spheres induced by theta polymers. Fur- 
thermore, since the depletion potentials are the domi- 
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FIG. 13: Depletion potentials, calculated via Eq. I|12^ . for 
p/p* = 0.26 and p/p* ~ 1.01. This simple theory compares 
well to direct simulations, taken from—, for good solvents. For 
theta polymers the depletion potential resembles that of ideal 
polymers: at p/p* = 0.26 they are virtually inistinguishable; 
at p/p* — 1.01 there are slight differences. 

nant determinant of depletion induced phase separation 
for Rg/Rc ^ l^^i^^ , we would expect the phase-behavior 
of theta polymers to closely follow that of ideal poly- 
mers, at least for these size ratios. If there are slight de- 
viations, then Fig. ^| suggests that the binodals would 
shift toward those of polymers in good solvent. Recent 
experiments^ of colloids mixed with polymers in a near 
theta solvent indeed show behavior that more closely re- 
sembles that of ideal polymers than that of polymers in 
good solvent. 



to be expected^JJ-. Working out the detailed crossover 
behavior would require many more simulations at differ- 
ent L. However, in practice, most experimentally investi- 
gated polymer solutions do not reach the scaling regime 
either. Moreover, we do observe significant qualitative 
differences between theta polymers, and polymers with 
weaker monomer- monomer attraction. 

In contrast to polymers in good solvent, solutions of 
polymers in theta solvent are quite well described by ideal 
polymer theories throughout the entire dilute regime. 
This works best for bulk properties like the EOS, and 
slightly less well for interfacial properties such as the 
density profiles, adsorptions, and surface tension. How- 
ever, as predicted by thcory^-^, in the semi-dilute regime 
theta polymers also begin to to exhibit important devi- 
ations from ideal-polymers. Our simulations show be- 
havior which is consistent with that predicted by scaling 
theories, but the concentration range we were able to in- 
vestigate is not wide enough to unambiguously confirm 
the expected scaling behavior. We also suggest that, for 
large enough concentrations, the values of various prop- 
erties, including the EOS, the adsorption, and the surface 
tension, will be larger for theta polymers than for poly- 
mers in good solvent at the same reduced concentration 
p/p*. 

The close agreement between theories for ideal poly- 
mers, and our simulations of theta polymers in the dilute 
regime, suggests that the effective depletion pair poten- 
tials, and associated phase behavior, should also resemble 
those of ideal polymers, at least in the so-called "colloid 
limit", where Rg/Rc ^ 1- 



VII. CONCLUSIONS 

We have performed extensive Monte Carlo simulations 
to investigate the bulk and interfacial properties of poly- 
mers in good solvent, and polymers in the theta regime. 
For infinite length L, there would be a sharp transition 
between the behavior of polymers in good solvent, and 
those in a theta solvent. The simulations were carried 
out for L — 500, so that continuous cross-over effects are 
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